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The Pulsar Wind Nebula (PWN) 3C 58 is energized by one of the highest spin-down power pulsars known (5% 
of Crab pulsar) and it has been compared to the Crab Nebula due to their morphological similarities. This 
object was detected by Fermi-LAT with a spectrum extending beyond 100 GeV. We analyzed 81 hours of 3C 
58 data taken with the MAGIC telescopes and we detected VHE gamma-ray emission for the first time at TeV 
energies with a significance of 5.7 sigma and an integral flux of 0.65% C.U. above 1 TeV. The differential energy 
spectrum between 400 GeV and 10 TeV is well described by a power-law function d^/dE = f 0 (E/lTeV)~ r 
with f 0 = (2.0 d= OAstat ± 0.6sys)10~ 13 cm~ 2 s~ 1 TeV~ 1 and F = 2.4 ± 0.2 sta d= 0.2 sys. This leads 3C 58 to be 
the least luminous PWN ever detected at VHE and the one with the lowest flux at VHE to date. According 
to time-dependent models in which electrons up-scatter photon fields, the best representation favors a distance 
to the PWN of 2 kpc and FIR comparable to CMB photon fields. If we consider an unexpectedly high FIR 
density, the data can also be reproduced by models assuming a 3.2 kpc distance. A low magnetic field, far from 
equipartition, is required to explain the VHE data. Hadronic contribution from the hosting supernova remnant 
(SNR) requires unrealistic energy budget given the density of the medium, disfavoring cosmic ray acceleration 
in the SNR as origin of the VHE gamma-ray emission. 


1. General description 

The supernova remnant 3C 58 (SNR G130.7+3.1) 
has a flat radio spectrum and is brightest near the cen¬ 
ter, therefore it was classified as a pulsar wind nebula 
[PWN; [l|. It is centered on PSR J0205+6449, a pul¬ 
sar discovered in 2002 with the Chandra X-ray obser¬ 
vatory Q. It is widely assumed that 3C 58 is located 
at a distance of 3.2 kpc 13]. but recent H I measure¬ 
ments suggest a distance of 2 kpc [4]. The age of the 
system is estimated to be ^ 2.5 kyr [5[ from the PWN 
evolution and energetics, however this is a matter of 
debate. The pulsar has one of the highest spin-down 
powers known (E = 2.7xl0 37 erg s _1 ). The PWN 
has a size of 9'x6' in radio, infrared (IR), and X-rays 
@, 7,0,0. Its luminosity is L 0 .5 - 10 keV = 2.4 x 10 34 
erg s -1 in the X-ray band, which is more than 3 or¬ 
ders of magnitude lower than that of the Crab nebula 
0- 3C 58 has been compared with the Crab because 
the jet-torus structure is similar [8]. Because of these 
morphological similarities with the Crab nebula and 
its high spin-down power (5% of Crab), 3C 58 has his¬ 
torically been considered one of the PWNe most likely 
to emit 7 rays. 

The pulsar J0205+6449 has a period P=65.68 ms, 
a spin-down rate P = 1.93 x 10 -13 s s -1 , and a char¬ 
acteristic age of 5.38 kyr [2j. It was discovered by the 
Fermi -LAT in pulsed 7 rays. The measured energy 
flux is F 7> o.iGeV=(5.4±0.2)xlO -11 erg cm -2 s -1 with 
a luminosity of T 7 >o .i GeV =(2.4±0.1)xl0 34 erg s" 1 , 
assuming a distance for the pulsar of 1.95 kpc [ui 


The spectrum is well described by a power-law with an 
exponential cutoff at E cuto ff=1.6 GeV 0. No pulsed 
emission was detected at energies above 10 GeV [13]. 
In the off-peak region, defined as the region between 
the two 7-ray pulsed peaks (off-peak phase interval 
0=0.64-0.99), the Fermi Collaboration reported the 
detection of emission from 3C 58 [l2|. The reported 
energy flux is (1.75T0.68) x 10 -11 erg cm -2 s -1 and the 
differential energy spectrum between 100 MeV and 
316 GeV is well described by a power-law with pho¬ 
ton index 7 = 1.61 ±0.21. No hint of spatial exten¬ 
sion was reported at those energies. The association 
of the high-energy unpulsed steady emission with the 
PWN is favored, although an hadronic origin related 
to the associated SNR can not be ruled out. 3C 58 
was tagged as a potential TeV 7-ray source by the 
Fermi Collaboration [13]. 

The PWN 3C 58 was previously observed in the 
VHE 7-ray range by several IACTs. The Whip¬ 
ple telescope reported an integral flux upper limit 
of 1.31 xlO -11 cm -2 s -1 ~ 19 % C.U. at an energy 
threshold of 500 GeV [14], and VERITAS established 
upper limits at the level of 2.3 % C.U. above an en¬ 
ergy of 300 GeV [l5|. MAGIC-I observed the source 
in 2005 and established integral upper limits above 
110 GeV at the level of 7.7xl0 -12 cm -2 s -1 (~4 % 
C.U.) [H| . The improved sensitivity of the MAGIC 
telescopes with respect to previous observations and 
the Fermi -LAT results motivated us to perform deep 
VHE observations of the source. 
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Figure 1: Distribution of squared angular distance, 0 2 , between the reconstructed arrival directions of gamma-ray 
candidate events and the position of PSR J0205+6449 (red points ). 


2. MAGIC observations and results 

MAGIC observed 3C 58 in the period between 4 Au¬ 
gust 2013 to 5 January 2014 for 99 hours, and after 
quality cuts, 81 hours of the data were used for the 
analysis. The data were analyzed using the MARS 
analysis framework The source was observed at 
zenith angles between 36° and 52°. The data were 
taken in wobble-mode [18] pointing at four different 
positions situated 0.4° away from the source to evalu¬ 
ate the background simultaneously with 3C 58 obser¬ 
vations. 

The applied cuts yield an energy threshold of 420 
GeV. The significance of the signal, calculated with 
the LiMa formula, is 5.7cr, which establishes 3C 58 as 
a 7 -ray source. The 0 2 distribution is shown in Figure 
[TJ As the five OFF positions were taken for each of 
the wobble positions, the OFF histograms were re¬ 
weighted depending on the time taken on each wobble 
position. 

We show in Figure [ 2 ] the relative flux (ex¬ 
cess/background) skymap, produced using the same 
cuts as for the 0 2 calculation. The TS significance, 
which is the LiMa significance applied on a smoothed 
and modeled background estimate, is higher than 6 
at the position of the pulsar PSR J0205+6449. The 
excess of the VHE skymap was fit with a Gaussian 
function. The best-fit position is RA(J 2000 )= 2 h 05m 
31(09) at at(ll)*y*s ; DEC(J2000)=64° 51'(l) at at(lW 
This position is statistically deviant by 2a from the 


position of the pulsar, but is compatible with it at la 
if systematic errors are taken into account. In the bot¬ 
tom left of the image we show the point spread func¬ 
tion (PSF) of the smeared map at the corresponding 
energies, which is the result of the sum in quadrature 
of the instrumental angular resolution and the applied 
smearing (4.7' radius, at the analysis energy thresh¬ 
old). The extension of the signal is compatible with 
the instrument PSF. The VLA contours are coincident 
with the detected 7 -ray excess. 

Figure [3] shows the energy spectrum for the MAGIC 
data, together with published predictions for the 
gamma-ray emission from several authors, and two 
spectra obtained with three years of Fermi -LAT data, 
which were retrieved from the Ferrm-LAT second 
pulsar-catalog [2PC, [ 12 | and the Fermi high-energy 
LAT catalog [1FHL, [PJ. The 1FHL catalog used 
events from the Pass 7 Clean class , which provides 
a substantial reduction of residual cosmic-ray back¬ 
ground above 10 GeV, at the expense of a slightly 
smaller collection area, compared with the^ Pass 7 
Source class that was adopted for 2PC 0. The 
two 7 -ray spectra from 3C58 reported in the 2PC 
and 1FHL catalogs agree within statistical uncer¬ 
tainties. The differential energy spectrum of the 
source is well fit by a single power-law function 
d(t>/dE=f 0 (E/l TeV)- r with f 0 = (2.0 ± 0.4 stat ± 
0.6 S y S )10 _ 13 cm _ 2 s _ 1 TeV _1 , T = 2.4 ± 0.2 stat ± 0.2 sys 
and x 2 =0.04/2. The systematic errors were estimated 
from the MAGIC performance paper j 2 l| including 
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Figure 2: Relative flux (excess/background) map for MAGIC observations. The cyan circle indicates the position of 
PSR J0205+6449 and the black cross shows the fitted centroid of the MAGIC image with its statistical uncertainty. In 
green we plot the contour levels for the TS starting at 4 and increasing in steps of 1. The magenta contours represent 
the VLA flux at 1.4 GHz 0 , starting at 0.25 Jy and increasing in steps of 0.25 Jy. 


the upgraded telescope performances. The integral 
flux above 1 TeV is Fe> i TeV = 1.4 x 10 - 13 cm - 2 s -1 . 
Taking into account a distance of 2 kpc, the lumi¬ 
nosity of the source above 1 TeV is L 1: e> i TeV = 
(3.0 zb 1.1) xlO 32 ^ erg s _1 , where d 2 is the distance 
normalized to 2 kpc. 


3. Discussion 

Several models have been proposed that predict the 
VHE 7 -ray emission of PWN 3C 58. 

[25] presented a one zone model of the spectral evo¬ 
lution of PWNe and applied it to 3C 58 using a dis¬ 
tance of 3.2 kpc. The VHE emission from this model 
consists of IC scattering of CMB photons and optical- 
to-IR photons, and also of pion decay. The flux of 7 
rays above 400 GeV predicted by this model is about 
an order of magnitude lower than the observation. 

[ 22 ] proposed a time-dependent model in which 
positrons gain energy in the process of resonant scat¬ 
tering by heavy nuclei. The VHE emission is pro¬ 
duced by IC scattering of leptons off CMB, IR, and 
synchrotron photons and by the decay of pions due 
to the interaction of nuclei with the matter of the 
nebula. The age of 3C 58 is assumed to be 5 kyr, 
using a distance of 3.2 kpc and an expansion ve¬ 


locity of 1000 km s -1 . According to this model, 
the predicted integral flux above 400 GeV is ^ 10“ 13 
cm - 2 s -1 , while the integral flux above 420 GeV mea¬ 
sured here is 5x 10 _ 13 cm _ 2 s _1 . Calculations by [26], 
using the same model with an initial expansion veloc¬ 
ity of 2000 km s -1 and considering IC scattering only 
from the CMB, are consistent with the observed spec¬ 
trum. However, the magnetic field derived in this case 
is B ^14/iG and it underestimates the radio emission 
of the nebula, although a more complex spectral shape 
might account for the radio nebula emission. 

[27 ] developed a time-dependent model of the spec¬ 
tral evolution of PWN including synchrotron emis¬ 
sion, synchrotron self-Compton, and IC. They evolved 
the electron energy distribution using an advective dif¬ 
ferential equation. To calculate the observability of 3C 
58 at TeV energies they assumed a distance of 2 kpc 
and two different ages: 2.5 kyr and 1 kyr [23j. For the 
2.5 kyr age, they obtained a magnetic field B ~ 17 /iG, 
while for an age of 1 kyr, the magnetic field obtained 
is B=40 /iG. The emission predicted by this model is 
closer to the Fermi result for an age of 2.5 kyr. 

[ 2 8| presented a different time-dependent leptonic 
diffusion-loss equation model without approxima¬ 
tions, including synchrotron emission, synchrotron 
self-Compton, IC, and bremsstrahlung. They as¬ 
sumed a distance of 3.2 kpc and an age of 2.5 kyr 
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Figure 3: 3C 58 spectral energy distribution in the range between 0.1 GeV and 20 TeV. Red circles are the VHE points 
reported in this work. The best-fit function is drawn in red and the systematic uncertainty is represented by the yellow 
shaded area. Black squares and black arrows are taken from the Fermi -LAT second pulsar-catalog results [12]. Blue 
squares are taken from the Fermi high-energy LAT catalog 0. The magenta line is the SED prediction for 3C 58 
taken from Figure 10 of [_22]. The clear green dashed-dotted line is the SED predicted by [33], assuming an age of 1 
kyr, and the dark green dotted line is the prediction from the same paper, assuming an age of 2.5 kyr. The blue dashed 
line represents the SED predicted by 0 assuming that the Galactic FIR background is high enough to reach a flux 
detectable by the MAGIC sensitivity in 50h. 


to calculate the observability of 3C 58 at high ener¬ 
gies [24]. The predicted emission, without considering 
any additional photon source other than the CMB, is 
more than an order of magnitude lower than the flux 
reported here. It predicts VHE emission detectable 
by MAGIC in 50 hours for an FIR-dominated photon 
background with an energy density of 5 eV/cm 3 . This 
would be more than one order of magnitude higher 
than the local IR density in the Galactic background 
radiation model used in GALPROP [~ 0.2 eV cm -3 ; 
[29]. The magnetic field derived from this model is 35 
fiG. To reproduce the observations, a large FIR back¬ 
ground or a revised distance to the PWN of 2 kpc are 
required. In the first case, a nearby star or the SNR it¬ 
self might provide the necessary FIR targets, although 
no detection of an enhancement has been found in the 
direction of the PWN. As we mentioned in Sec. [lj a 
distance of 2 kpc has recently been proposed by [4] 
based on the recent H I measurements of the Cana¬ 
dian Galactic Plane Survey. At this distance, a lower 
photon density is required to fit the VHE data. 

We have shown different time-dependent models in 
this section that predict the VHE emission of 3C 58. 
The SED predicted by them are shown in Figure [3j 


They use different assumptions for the evolution of 
the PWN and its emission. [25] divided the evolution 
of the SNR into phases and modeled the PWN evo¬ 
lution inside it. In [22|] model, nuclei play an impor¬ 
tant role in accelerating particles inside the PWN. [24| 
and J27] modeled the evolution of the particle distri¬ 
bution by solving the diffusion-loss equation. [24] fully 
solved the diffusion-loss equation, while [ 27 J neglected 
an escape term in the equation as an approximation. 
Another difference between these latter two models is 
that [27] took synchrotron emission, synchrotron self- 
Compton and IC into account, while (24[ also consider 
the bremsstrahlung. The models that fit the 7 -ray 
data derived a low magnetic field, far from equiparti- 
tion, very low for a young PWN, but comparable with 
the value derived by [9] using other data. 


4. Conclusions 

We have for the first time detected VHE 7 rays up 
to TeV energies from the PWN 3C 58. Following the 
assumptions in [30j , it is highly unlikely that the mea¬ 
sured flux comes from hadronic emission of the SNR. 
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The measured luminosity and flux make 3C 58 into 
an exceptional object. It is the weakest VHE PWN 
detected to date, a fact that attests to the sensitiv¬ 
ity of MAGIC. On the other hand, it is also the least 
luminous VHE PWN, far less luminous than the orig¬ 
inal expectations. Its ration Lvhe/F ~ 10 -5 is the 
lowest measured, similar to Crab, which makes into a 
very inefficient 7 -ray emitter. Only a closer distance 
of 2 kpc or a high local FIR photon density can qual¬ 
itatively reproduce the multi wavelength data of this 
object in the published models. Since the high FIR 
density is unexpected, the closer distance with FIR 
photon density comparable with the averaged value 
in the Galaxy is favored. The models that fit the 7 - 
ray data derived magnetic fields which are very far 
from equipartition. 
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